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Abstract 



m 
(N 

J> ■ We present an exact expression for the evolution of the wavefunction of N 

o 

*^ ■ interacting atoms in an arbitrarily time-dependent, d-dimensional parabolic 

■ trap potential uj{t). The interaction potential between atoms is taken to be of 

0^ . the form .^/r^ with ^ > 0. For a constant trap potential uj{t) = ljq, we find an 



exact, infinite set of relative mode excitations. These excitations are relevant 



B 

. to the linear response of the system; they are universal in that their frequen- 

O ■ cies are independent of the initial state of the system (e.g. Bose-Einstein 



condensate), the strength of the atom-atom interaction, the dimensionality 
d of the trap and the number of atoms A^. The time evolution of the sys- 
tem for general uj{t) derives entirely from the solution to the corresponding 
classical ID single-particle problem. An analytic expression for the frequency 
response of the A^-atom cluster is given in terms of uj{t). We consider the im- 
portant example of a sinusoidally-varying trap perturbation. Our treatment, 
being exact, spans the 'linear' and 'non-linear' regimes. Certain features of 
the response spectrum are found to be insensitive to interaction strength and 
atom number. 

PACS numbers: 32.80.Pj, 42.50. Vk,67.65.+z 
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I. INTRODUCTION 

Atomic traps have been the subject of much recent research activity, both theoretical |1| 
and experimental 0. The finite atom number and reduced effective dimensionality 
(i < 3 in atomic traps opens up a fascinating research area of confined A^-atom states. A 
singularly important milestone in this field has been the recent experimental observation of 
a Bose-Einstein condensate in an effectively two-dimensional {d = 2) trap 0. It is to be 
expected that the application of time-dependent perturbations to such traps will provide 
useful information as to the properties of the A^-atom system. 

According to the preparation of the atomic trap, the confinement length scales in the 
three spatial directions {L^, Ly and L^) can, in principle, be quite different yielding highly 
anisotropic traps. Within a simple particle-in-a-box picture the single-particle energy level 
spacing AE ~ L^^', it follows that for >> Ly, the atoms will be stuck in the lowest y, z 
subbands hence freezing out the y, z degrees of freedom. The anisotropic three-dimensional 
trap is now effectively quasi-one-dimensional in that the atoms only have significant free- 
dom along the x-direction. Similarly if ~ -^j/ >> -^z the trap is effectively quasi-two- 
dimensional; if Ly Lz the trap becomes quasi-three-dimensional. 

A variety of theoretical predictions have been made concerning the properties of atomic 
traps. For example, Burnett and co-workers have obtained very interesting results concerning 
the linear and non-linear responses of trapped neutral atoms in a Bose-Einstein condensate 
by solving a time-dependent Gross-Pitaevskii equation for the condensate order-parameter 
or 'wavefunction' [^. Important open questions concern the precise dependence of the A^- 
atom quantum state on the number of atoms A^ [^], and on the strength of the interaction 
between atoms 0. 

Given the general intractability of the A^-body problem, most theoretical approaches have 
been numerical and tend to be based on macroscopic, phenomenological models or mean- 
field approximations. Our results in this paper complement such studies. In particular 
we provide analytic results for the response of an arbitrary number A^ of trapped, neutral 
atoms by studying a microscopic model Hamiltonian. We provide closed-form expressions 
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for the trap properties as a function of the trap parameters, thereby yielding insight into the 
competition between the confinement and atom-atom interactions. The analytic tractability 
of our model is made possible through a combination of the parabolic form for the confining 
potential in the magnetic traps of interest [Q], and an inverse-square interaction between 
atoms, i.e. l/r" with n = 2. The actual short-range repulsive interaction between neutral 
atoms is likely to be better fitted by a larger n; however, the general features of our results 
should be qualitatively similar for other short-range interactions involving n > 2. 

The paper is organized as follows. In Sec. II we discuss the excitation spectrum obtained 
from the time-independent A^-atom Schrodinger equation for an isolated trap. These excita- 
tions are relevant for describing the response of the trap to a weak, sinusoidal perturbation 
(i.e. linear response) . In Sec. Ill we consider a class of time-dependent perturbations having 
arbitrary strength and form, and derive the exact, time-dependent A^-atom wavefunction. 
As an example, the exact, non-linear time-dependent response to a sinusoidal perturbation 
is developed in Sec. IV; this response is shown to be largely determined by solutions of 
the well-known Mathieu equation. Non-linear features in the system response are identified 
and explained. Throughout the paper we will focus on quasi-two-dimensional traps, since 
these systems are attracting much experimental interest. The generalization of these results 
to quasi-one- dimensional and isotropic, three-dimensional traps is straightforward and the 
corresponding results are discussed briefiy. We will also focus on the case of repulsive in- 
teractions between atoms. The extent to which our results are also applicable to attractive 
atom-atom interactions is discussed in Sec. V, together with the conclusions. We will com- 
ment throughout on the extent to which our analytic results can reproduce similar physics 
to recent numerical work. 

II. ELEMENTARY EXCITATIONS OF AN ISOLATED TRAP 

In this section we focus on the elementary excitations of an isolated trap. Such excitations 
govern the linear response of the system in the presence of a weak, sinusoidal perturbation. 
We will specifically obtain an exact, infinite set of relative mode excitations for the A^-atom 
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system. 

The governing microscopic equation for the energy spectrum of the atom trap is the N- 
atom non-relativistic Schrodinger equation with a d-dimensional confining potential Vc{r); 
this is given by if^ = with 

^ = E [^P' + M^i)] + E ^/(r. - r,) , (1) 

i i<j 

where the momentum and position associated with the i'th atom are given by pj and r^. The 
magnetic traps of interest have approximately parabolic confining potentials Vc- We will 
assume that the trap potential for a d-dimensional trap is a d-dimensional isotropic harmonic 
oscillator, ie. Vc — ^mwlrf. We also assume the atom-atom interaction is translationally- 
invariant. The exact functional form of the two-body interaction Vi{Yi — Yj) between neutral 
atoms is unknown; here we focus on a net repulsive interaction between atoms of the form 
^(|rj — rj|)~" where n — 2 and ^ > 0. Note that for n > 2, the interaction is sufficiently 
'short-range' that we would expect similar behavior to emerge for a range of n. At the end of 
the paper we briefly discuss how our results carry over to an attractive interaction between 
atoms (i.e. ^ < 0). 

We employ standard Jacobi coordinates Xj (i = 0, 1, . . . , A?" — 1) where Xq = ;^ J2j 
(center-of-mass), Xi = ■y|(r2 — ri), X2 = — r^) etc. together with their conjugate 

momenta Pj. The center-of-mass motion decouples, H — HcmO^o) + -f^rei({Xi>o}), hence 
E — Eqm + -E'rei and ^ = V'civiV'rei- The exact eigenstates V'cm of Hqm and eigenenergies 
Eqm are identical to those of a single particle in a parabolic potential and are well-known. 
We will concentrate on the relative motion since it is here that effects of atom-atom correla- 
tions manifest themselves. The non-trivial problem is to solve the relative motion equation 
ifj-ei'^rei = -^'rei'^/'rei- We transform the relative coordinates {Xj>o} to standard hyperspher- 
ical coordinates: Xj = r^Ufs^"^ cosaj+i)smaie^^' with r > and < ctj < | (cti = |). 
Physically, the hyperradius r is related to the root-mean-square atom-atom separation: 

Af-l 1 N 

r'-E^^-i^ E (r.-r.f (2) 

j=l (i>i') i=2 
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As mentioned above, we will focus on d = 2. The exact eigenstates of i^rei have the form 
ipj-ei = R{r)F{Cl) where Cl denotes the (2A^ — 3) remaining {6, a} variables; R{r) and F{fl) 
are solutions of the hyperradial and hyperangular equations respectively. The hyperradial 
equation is 

d' 2N-3d 7(7 + 2iV - 4) r\2mE^ 

+ ^-Tr 7- I + ^^^^^^^ = ° 

where /q = h{muQ)~^; the parameter 7 is related to the eigenvalue of the tuo-independent 
hyperangular equation (see below). Equation (3) can be solved exactly yielding 

E,ei = hcuoi2n + -f + N-l) (4) 

where n is any positive integer or zero and 

^n(r) = (f rLr^-^(ie"4 . (5) 

where L is the Laguerre Polynomial. Equation (4) provides an exact (and infinite) set of 
relative mode excitations at frequencies 2ujQAn, irrespective of the initial eigenstate of the 
A^-atom system (e.g. Bose-Einstein condensate). These are 'breathing' modes; we note that 
they are very similar in frequency to the set of modes found to occur near multiples of 2uo 
in recent numerical calculations based on the Gross-Pitaevskii equation [Q. 
The cUQ-independent hyperangular equation determining 7 is given by 

[0^ + ^-^VmF{n) = [7(7 + 2N- A)]F{n) (6) 

where 

d' [2Ar - 6 - (2iV - 4)cos2a^] d , , d' 

'^N = -Tr^^ ^ 7^ h sec a7ve^„i - cosec aAr--2- . (7) 

oaff sm2aAr daN o9% 

The quantity V{VL) represents the hyperangular part of the atom-atom interaction and is 
given by 

m=-'ET-^. (8) 



Since the atom-atom interaction only depends on absolute relative coordinates, V{fl) com- 
mutes with the relative angular momentum. The relative angular momentum J therefore 
remains a good quantum number, hence we can introduce a further Jacobi transformation 
of the relative motion angles {Oi}: in particular 9' — Sili^ ^i, ^ — ^(^^ ~ ^2) etc. The 
hyperangular equation now depends on just {2N — A) remaining hyperangles, i.e. the hypcr- 
angles fl excluding 9'. We denote these (2A^ — 4) hyperangles as fl. The exact eigenstates 
of //rei now have the form V'rei = e'-^^' R(r)G{n) (i.e. F(Q) = e^-^^'G'(Q)). 

Unfortunately the hyperangular equation (Eq. (6)) does not admit complete exact so- 
lutions for 7 and F{fl), or equivalently G{fl). However some insight into the properties 
of these eigenstates can be gained as follows. Without loss of generahty, we can choose to 
rewrite the eigenvalue of the hyperangular equation (Eq. (6)) in terms of a new quantity e 
as 

6=^[7(7 + 27V-4)-(^)^-J^] (9) 

where Kiass is the potential energy of the minimum-energy configuration for classical atoms 
in the trap (i.e. the Wigner solid). This quantity e is useful in that it isolates purely 
quantum-mechanical contributions to 7, i.e. it does not contain either the classical potential 
energy of the atoms or their rotational kinetic energy. It is straightforward to show that 
Kiass ^^ujo and hence e (like 7) is independent of ujq. We can now recast the exact relative 
energy expression as 

£;,ei = nu;o[2n + {[N - 2]' + j' + + S)^ + 1] ■ (10) 

ncuo h 

Erei only depends on particle statistics through e. As ^ — > 0, Ej-^i — > Kiass and e — > 
0. Physically, e includes the zero-point energy in fi-space associated with the quantum- 
mechanical spread of G(0) about the hyperangles fl corresponding to the classical, minimum 
energy configuration of the A'"-atom system (Wigner solid). In the limit of fairly weak 
interatomic interactions (i.e. small ^), the spread in G{Q) and hence the magnitude of the 
zero-point energy e will be large. For any given ^, the spread in G(0) and hence magnitude of 
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e will also depend on total wavefunction symmetry requirements; e will in general be smaller 
for a bosonic ground state than for a fermionic ground state due to the lower number of 
nodes in G{Q) for bosons. 

The above analysis which yielded the exact excitation frequencies 2ujQAn for quasi-two- 
dimensional {d = 2) traps carries over to quasi-one-dimensional traps {d = 1) and isotropic 
three-dimensional traps {d = 3). The same breathing- mode frequencies 2uJoAn are obtained 
in each case. For d = 1 the quantity 7 in Eq. (4), and hence the entire energy spectrum Ej-^i, 
can actually be obtained analytically. This model was first solved exactly by Calogero 0. 
The exact A^-atom energy levels for a quasi-one-dimensional trap are given hj E = Ecu+Ej-ei 
with 

i?rci = huJoiliN - 1) + liV(iV - l)(r + i) + k) (11) 

where r= + and is a positive integer (N.B. k ^ 1). The hyperradial breathing 

modes 2huJoAn are included in this relative mode spectrum via the integer k values. We 
note that Eq. (11), while exact for non-zero ^, does not yield the full energy level spectrum 
for non-interacting atoms in the limit — 0. This feature, which is discussed explicitly 
in Ref. [7], results from the singular nature of the inverse- square interaction when applied 
in one dimension - it does not arise in two or three dimensions. 



III. DYNAMICS OF A PERTURBED TRAP 

In this section we extend the above results to include a time-dependent trap potential 

( for t < 

u\t) = (12) 
[ f{t) for t > 

The treatment will be exact, analytic and applicable to any function f{t) with arbitrarily 
large magnitude, i.e. we implicitly include all non-linear effects to all orders. Specifically, we 
will provide the time-dependent versions of the hyperradial breathing-mode wavefunctions 
given in Eq. (5). These will be found to depend only on the solution of a classical one- 
dimensional oscillator with trap potential uj'^{t). Again we will specifically consider a quasi- 



two-dimensional {d = 2) trap, but the method is equally applicable to both d = 1 and 
d = 3. 

The time-dependent Schrodinger equation is given by H'^ = ih{d'^ /dt). In the presence 
of the time-dependent trap potential discussed above, the separation of the center-of-mass 
and relative motion is still exact, \E' = ipcM{t)'4've\{t)- It was found in Sec. II that the 
hyperangular equation is independent of tuo, hence those parts of the relative wavefunction 
iprei that derive from this equation remain time-independent. The time-dependence of the 
total wavefunction \E' is therefore only contained in the center-of-mass and hyperradial parts. 
Here we consider explicitly ipreiit) and hence the hyperradial part R{r,t). The solution for 
the trivial center-of-mass part ipcuit) is exactly analogous. Following we construct the 
generating function 

oo 

g{z,r,t) = Y,R^{r,t)z^ (13) 

n=0 

where Rn{r,t) are the solutions of the time-dependent hyperradial equation. Because the 
trap potential is constant for t < 0, we can employ Eq. (5) to obtain 

oo 2 _ri 

g{z,r,t < 0) = E (^)"^r''~^(72)e "^^^ (14) 

n=0 '0 '0 

This can be written in closed form using the identity 

oo 

$:^"L:i(y)^(i-^r('^+v^) (15) 

n=0 

(see Ref. ^) so that 

g{z,r,t< 0) = (f )^e^^^^'(l -2)-("+^) (16) 

where a = 7 + — 2. We now make the ansatz 

^?(^,r,t>0) = a(^,t)(f)^e"'(^'*)^' (17) 

to 

which can be shown to satisfy the time- dependent hyperradial equation 

, d'^ 2N-3d 7(7 + 2A^-4) ,m,^ 2/ n 2x / x 2im d , 



and the t — boundary condition in Eq. (16), provided 

a{z,t) = [rj{t)]-^''+^^exp[2ie{t){a + - z exp[2i^(i)])-("+^) (19) 

and 

"'^^'^^ = " ^'^f^^^^ " ' ea;p[2.^(t)])-^) (20) 

where r]{t) = \ri{t)\e^^^^^ solves the classical one- dimensional oscillator 

rj'it) + f{t)rj{t) = (21) 

with boundary conditions 77(0) = 1 and 77(0) = —icuo- We may then expand g{z,r,t) using 
the identity (15), and compare coefficients of with the defining Eq. (13) to obtain the 
desired time-dependent wavefunctions (unnormalised) : 

R^{r,t > 0) = |r7(t)|i-Vexp[i(^(i)(2n + a + 1) + £-±\r^{t)\')]e-b' L^iy') . (22) 

where y = In the static case this expression reduces to equation (5). If for a particular 
f{t) we are able to solve Eq. (21), then Eq. (22) provides a complete description for 
the evolution of the initially stationary hyperradial state as the trap becomes time- 
dependent for i > 0. As noted previously, since the remaining hyperangular part of V'rei is 
time-independent, Eq. (22) together with the (exactly analogous) expression for the center- 
of-mass entirely determines the time-evolution of the total wavefunction ^ = '0CM(^)'0rei(^)- 
We now comment on the significance of our results, particularly with respect to the Bose- 
condensed atomic gas. Our results for the breathing modes in Sec. II and for the dynamical 
response in Sec. Ill are valid for all strengths of the atom-atom interaction (^), all trap sizes 
(a;o) and all numbers of atoms (N). They are independent of the precise form of the initial 
state of the system and are therefore true in both the Bose-condensed and non-condensed 
regimes. Since we do not know 7 (or F{fl)) explicitly, our model does not allow us to predict 
a temperature at which the gas will be substantially Bose-condensed. This is consistent with 
the expectation that the condensation temperature has a complicated dependence on a;o, C 
and N. 
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In deriving the time-dependent wavefunction of Sec. Ill, we took the A^-atom system to 
be in a single eigenstate of the time-independent Schrodinger equation up until the moment 
the perturbation is turned on (i.e. for t < 0). The linearity of the Schrodinger equation 
implies that any initial wavefunction consisting of a sum of eigenstates will evolve as the 
sum of its parts, and may therefore be written as a sum of terms such as Eq. (22). In Sec. 
IV below, the calculated response spectrum is also derived under the assumption of a single 
eigenstate for t < 0, however the initial state can be generalized to a sum of eigenstates 
without altering the response spectrum. The results we present in this paper are therefore 
remarkably general in their applicability. 



IV. EXAMPLE: SINUSOIDAL PERTURBATION 

As an illustration of the applicability of our formalism, we now take a particular f{t) 
and employ Eq. (22) to calculate the frequency response spectrum of the system. We find 
that the spectrum exhibits some features which are essentially independent of 7. Since the 
^-dependence only enters the response through 7, it follows that any observables found to 
be substantially independent of 7 are generic to all atom traps independent of the strength 
of the atom- atom interaction ^. 

We consider the specific case of a sinusoidal perturbation turned on at t = 0: 

( for t < 

u;\t) = (23) 
[/(t) =c<;2-cj2(l-cos(2fit)) fort>0 

A sketch of this function appears at the top of Fig. 1. It would be quite possible to produce 

this trap perturbation experimentally. The solution to the classical equation (Eq. (21)) is a 

general Mathieu function, which may be written in the form 

00 00 
ri(t) = Ae"^' C2ne^'''^' + ^e"^™ ^ CsnC-^^'^™ (24) 

n=— 00 n=— 00 

where and {c2n} are determined by a set of simultaneous equations (see Ref. 0). 

In the main part of Fig. 1 we show the behavior of /i as a function of the parameters uq, 
Ui and Q. This behavior is non-trivial; in certain regions (shown white) of the parameter 
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space /i is purely imaginary while in others (shown dark) it has a real part. We define 
(3 = Jm{yu}. A non-zero real part of indicates that the classical particle is resonating 
with the oscillating trap; the particle's oscillations then become infinitely large as t ^ oo. 
The corresponding effect on the quantum mechanical system, which depends on time only 
through ?7(t), will be an increase in energy and a decrease in localization of i?(r, t) and hence 
V'rei- This spreading in ■j/'rei implies an increase in the average atom-atom separation, and 
will lead to atoms escaping from any realistic trap having a finite depth. For small cui, ix 
has a real component only when the perturbing frequency Vt is equal to an integer fraction 
^ of the trap potential ujq. We emphasize that our treatment is exact for any amplitude ujI 
of the sinusoidal perturbation. Indeed the theory remains sound even when ojI > ^uJq, in 
which case the trap becomes repulsive for part of each oscillation. 

In the remainder of this section we concern ourselves with solutions which are 'stable', i.e. 
points in the parameter space corresponding to a purely imaginary /i. As in the numerical 
studies of Ref. Q , we will calculate the frequency response spectrum of the atomic gas by 
Fourier analyzing the temporal variation of the gas density, the 'single-particle-density', at 
a fixed point in space. We define the single-particle-density at a point S as 



where \& is the total wavefunction at time t which, in our case, exactly separates into center- 
of-mass, hyperradial and angular parts. If we set S = (i.e. the trap center), we find 



The superscript on ipcuit) indicates that under the S = constraint, the center-of-mass 
coordinate is no longer independent of the other coordinates. Examining the form of ipcM 
and R we find that the integral in Eq. (26) may be rewritten in time-independent form 
using the scaled variable y = |^^[^|^^ , as employed in Eq. (22). This transformation to a 
time-independent integral results in an external factor |?7(t)|~^, thus the density is found to 
have a simple time-dependence: 




(25) 




(26) 
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d(S = 0,t) oc |r/(t)|-2 . (27) 

Interestingly, this expression (Eq. (27)) is completely general for any time- varying trap 
co'^(t); here we have chosen to consider a sinusoidal form for which the corresponding 77 (i) 
may be found analytically. The expression for (i(S = 0, t) is independent of both and the 
unknown quantity 7, thus the peak positions in the Fourier transform must be independent 
of variations either in the atom-atom interaction strength or in the number of atoms 
confined in the trap. Figure 2 displays the Fourier transform of Eq. (27) calculated with a 
particular choice of the parameters ojq and ui. We see a number of sharply defined peaks; 
the frequencies at which these peaks occur are all of the form nVt + mVL(3 where n and m 
are integers. The value of (3 may be found from Fig. 1; for the present choice of parameters 
(3 = 1.409. 

Figure 1 also shows that when the perturbation amplitude ui is small then ~ ^, where 
u'^ = Uq — ujf is the mean trap confinement. Thus for a weak-to-moderate perturbation, the 
peaks in the response lie at sums and differences of u and the driving frequency Q; i.e. at 
nQ + muj. Figure 1 shows us exactly how (3 (and hence the spectrum) changes as we move 
to the strongly non-linear regime {ui ^ u) . The parameters chosen for Fig. 2 correspond to 
a 'moderate' perturbation {wf = ClSScUg); the peaks marked n/3 in Fig. 2 deviate from nujQ 
by about 8%, and from nui by about 0.3% . Although a direct comparison with the non- 
linear response of Ref. is not practical because of the different form of the perturbation 
(i.e. white noise) it is interesting to note that the same types of non-linear phenomena are 
observed; in particular, harmonic generation and frequency mixing. As noted in Ref. [Q, 
such non-linear effects open up the possibility of non-linear atom optics based on coherent 
matter waves. 

V. CONCLUSION 

We have presented an analytically-solvable model of the quantum-mechanical time- 
evolution of A^ interacting atoms in an arbitrarily time-dependent, d-dimensional parabolic 



12 



trap potential u!{t). The solution allows us to determine the resonance frequencies of the 
system in response to both weak (linear response) and strong (non-linear response) pertur- 
bations. An exact, infinite set of relative mode excitations were found for a constant trap 
potential u!{t) — cuq; these excitations are universal in that their frequencies are independent 
of the initial state of the system (e.g. Bose- Einstein condensate), the strength ^ of the atom- 
atom interaction, the dimensionality d of the trap and the number of atoms N. The specific 
example of a sinusoidally- varying trap perturbation was employed to demonstrate the for- 
malism. Certain features in the system response were found to be insensitive to interaction 
strength and atom number. 

Finally we will discuss the extent to which our model can describe an attractive interac- 
tion between atoms (i.e. ^ < 0). The general formalism using hyperspherical coordinates is 
still valid; the difference is that e (see Eq. (9)) can now be negative, hence there is a possi- 
bility of collapse of the A^-atom system if the attraction is too strong (i.e. the atoms collapse 
to form an infinitely dense gas) . Such a complete collapse is an artefact of the inverse-square 
attractive interaction. It could be prevented by inserting an additional hard-core repulsion 
into the model, but the solvabihty of the model would be lost. Consider first a pair of atoms 
{N — 2) with an attractive interaction (^ < 0). In a quasi-one-dimensional trap {d — 1), 
it can easily be shown that two-atom collapse occurs if the strength of the attractive inter- 
action 1^1 > only weak attractive interactions where |C| < ^ therefore described 
by the present theory. The reason that a pair of atoms with a weak attraction can sup- 
port a non-collapsed ground state in the trap is that the negative potential energy is offset 
by a large, positive kinetic energy as a result of the uncertainty principle. In an isotropic 
three-dimensional trap {d — 3), two-atom collapse will occur if the strength of the attractive 
interaction |^| > ^(1 + 4/(/ -|- 1)) where / is the relative angular momentum of the two 
atoms. Hence for I — 0, the same condition holds as for d = 1. For quasi- two-dimensional 
traps, it turns out that two-atom collapse will occur for |^| > where is the component 
of relative angular momentum perpendicular to the two-dimensional plane. Therefore, for 
Iz — an arbitrarily small attractive potential will cause two-atom collapse. For atoms 
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in a quasi-one- dimensional trap, the same condition holds as for two atoms. For quasi-two 
and three-dimensional traps, the exact condition is not known since the problem is not com- 
pletely solvable (7 and hence e are not known exactly for > 2). It remains to be seen, 
therefore, whether the increase in kinetic energy during collapse is sufficient to offset the 
decrease in potential energy. If such a condition were to be satisfied, the A^-atom system 
with attractive interactions would form a ground state that was qualitatively different from 
the A^-atom system with repulsive interactions. As a matter of interest, we note that a 
suggestion of possible alternative states for A^ atoms with attractive interactions (^Li) has 
recently appeared |1T0| . 
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useful discussions. S.C.B. is supported by an EPSRC studentship. Partial funding was also 
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Figure Captions 



Figure 1. Top: The trap potential u;'^{t) (cf. Eq.(23)). Bottom: Contour plot showing the 
behavior of the constant /j, in the Mathieu function ri{t) (Eq. (21)) as a function of the 
trap potential parameters. This quantity /j, is complex, /j, = Re{ii} + ij3. White regions: 
Re{ii\ — and the contours are lines of constant /3; for a;i — > 0, /9 — > ^. Dark regions: 
i?e{/i} > and the contours are lines of constant Re{ii) in increments of 0.2. 

Figure 2. Frequency response spectrum showing Fourier amplitude as a function of frequency 
response for a trap with parameters uj\ = 0.138a;o and fl^ = 0.431a;Q. For these parameters 
(3 — 1.41; the (3 for any particular trap parameters (a;o,a;i,Jl) may be found using Fig. 1. 
Note that the two central lines (shown as cut-off) exceed the height of the graph by a factor 
of 3. 
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